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PROCEDE DE DETERAAINATION D'UNE STRUCTURE EN 3D A PARTIR 
D'UNE IMAGE EN 2D, NOTAMMENT LA STRUCTURE D'UN OS 

La prescntc invention conccrnc un procede de determination d'une structure en 
3D a partir d'une image en 2D, notamment celle d'un os. 

L'osteoporose est une maladie qui affecte en grand nombre surtout les femmes 
apres la menopause vers Tage de 50 ans sachant que cette maladie peut affecter 
toute personne et a tout age. Elle est caracterisee par une faible masse osseuse 
et une deterioration du tissu osseux. Ces degradations conduisent a des risques 
importants de fractures particulierement de la hanche, de la colonne vertebrale 
et du poignet. 

II existe bien sur un certain nombre de facteurs de risques mais ceci ne suf f it 

pas pour determiner la pro babi lite pour une personne de se f racturer un element 

de son squelette sachant qu'il n'y a aucun symptome prealable. 

Cette maladie est courante dans la population et devrait aller en s'etendont du 

fait de la sedentarisation notamment et du vieillissement de la population. 

Or le traitement des personnes af f ectees apres fracture est tres long et tres 

couteux car ce sont des so ins a long terme. Les consequences sont sou vent 

graves car conduisant a des invalidites voire a la mort des patients. 

Aussi est-il necessaire de prevenir ce type d' affection mais a condition de 

determiner le risque de fracture, ce qui est du seul ressort du praticien. 

Si Ton arrive a predire les risques, le praticien dispose alors de traitements par 

hormonotherapie et par prise d'autres produits de la pharmacopee comme la 

calcitonine ou les biphosphonates, en plus de conseils d' hygiene de vie pour 

tenter d'enrayer cette affection. 
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II existe dcs apparcillages pour determiner la densite osseuse appeles 
densitometres. 

On peut penser que des lors que Ton connaTt la densite osseuse, en comparant 
avec des courbes etalon, le praticien peut alors determiner les risques et etablir 
5 un diagnostic adapte et prescrire un traitement satisfaisant. 
C*est ce qui est couramment utilise. 

Par la demande de brevet US 5 774 520, on sait que la probabilite d'une 
fracture est directement liee a la densite osseuse. 

Les densitometres disponibles dans le commerce utilisent les rayons X ou les 
10 rayons gamma de faible energie. L'absorption etant proportionnelle a la densite 
de I'os, on peut done penser a correler les deux pour avoir un resultat 
satisfaisant permettant des comparaisons. Ce brevet prevoit un travail avec 
deux sources d* emission de photons en sorte de traiter la dispersion engendree 
par le fait que les matieres osseuses peuvent etre considerees comme 
15 homogenes mais que les autres matieres liquides, les tissus musculaires, les 
cartilages ne peuvent etre consideres comme une matiere unique, homogene. 
Un autre brevet US 6 385 283 utilise la densite mais combine cette mesure avec 
des images permettant de determiner les risques de fracture. Ces images sont 
realisees sur la colonne vertebrale du patient et Toperateur determine la 
20 presence d' amorces de fractures. De plus, comme dans Tart anterieur connu, 
ces informations sont comparees avec des donnees acquises au fur et a mesure 
sur de tres nombreux cas. Ce diagnostic est affine par le praticien a partir des 
antecedents du patient et par introduction des facteurs de risques 
supp lementai res . 

25 La demande de brevet WO 86/07531 propose de realiser une image d'un os 
donne, par exemple le calcaneum, a un age auquel le patient presente 
normalement sa plenitude de masse osseuse. Ensuite, le procede consiste a 
realiser des images successives de ce meme os pour permettre au praticien de 
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realiser des comparaisons ct de determiner revolution afin de conclure sur le 
degre de risque de fracture. 

Toutes ces methodes et autres dispositif s sont bases sur la seule mesure de la 
densite de la matiere osseuse. 
5 Or pour que le praticien arrive a un diagnostic satisfaisant, il faut qu'il ait des 
informations autres que la seule densite. 

La definition de I'osteoporose est donnee par T Organisation Mondiale de la 
Sante : ^'desordre caracterise par une faible masse osseuse et des alterations 
de la microarchitecture du tissu osseux, conduisant a augmenter la fragilite de 

10 I'os et en consequence a augmenter le risque de fracture " 

C*est ainsi qu'il existe des analyses de matiere osseuse in vivo plus sophistiquees 
a partir d' analyses de sang et d' urines faisant intervenir des marqueurs osseux 
mais ces analyses sont couteuses, contraignantes pour le patient et, de toutes 
les fagons, donnent des indications de composition de la matiere mais pas 

15 d ' indications sur sa structure. 

Or, on peut comprendre que la matiere osseuse presente une resistance 
differente suivant qu'il y a plus ou moins de matiere en quantite ce qui resulte 
de la mesure de densite. Par centre, pour une meme densite osseuse done pour 
une meme quantite de matiere, suivant T architecture de la structure osseuse, la 

20 resistance mecanique peut varier grandement. 

Ce qui est utile pour un praticien, ce sont les proprietes mecaniques de Tos d'un 
patient induisant directement les risques de fracture. Quand il y a diminution de 
la resistance osseuse, il y a conjugaison de deux facteurs, la baisse de la densite 
minerale osseuse DMO et la desalteration de la structure osseuse. 

25 A ce jour, il n' existe pas de moyens permettant de determiner la structure 
osseuse. On peut utiliser de T imager ie par resonance magnetique ou le scanner 
mais ce sont des examens lourds et couteux, surtout s'il faut effectuer un suivi 
longitudinal des patients. 
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De plus, la capture d' images en 3D a partir de moyens d*imagerle est pour 
rinstcnt limitee par les capacites de ces moyens car la resolution spatiale est 
superieure aux dimensions de la plupart des travees de liaison permettant de 
realiser la connectivite, c'est-a-dire les liaisons entre les noeuds du reseau en 
5 trois dimensions. 

Le procede selon la presente invention est maintenant decrit en detail, af in de 
permettre par Tanalyse d'une image 2D d'estimer les parametres mecaniques. 
Selon I'invention, le procede de determination de la resistance mecanique d'un os 
selon r invention, a partir d*une image 2D numerisee, obtenue par imagerie, se 

10 caracterise en ce que Ton realise une correlation entre la densite minerale 
osseuse DMO determinee a partir de cette image 2D par tout moyen adapte a ce 
type d'image et un parametre de structure obtenu a partir de la meme image 2D. 
Plus particulierement, on recourt a une correlation de type exponentielle. 
On utilise cette correlation associant la densite minerale osseuse et ledit 

15 parametre de structure pour determiner la contrainte ultime Cu de Tos. 

De f a^on plus precise, on determine le parametre de structure a , obtenu par la 
succession des etapes suivantes : 

a) choisir au hasard un pixel de T image 2D qui est a un niveau de gris h(0), 

b) choisir une droite issue de ce point ayant une direction egalement determinee 
20 au hasard, 

c) se deplacer d'une distance r sur cette droite, h(r) etant le niveau de gris de 
ce nouveau point, 

d) calculer la variance des niveaux de gris avec la formule : V(r) = [h(r) - h(0)]^, 

e) tracer la courbe associee a V(r) en echelle log-log, et 

25 f ) determiner la pente de cette courbe log-log qui represente ledit parametre a. 
Afin d'ameliorer la precision, on repete les etapes a) a d) un nombre de fois 
suff isamment grand pour faire converger la fonction variance V(r), moyennee 
sur r ensemble des iterations. 
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Selon une autre caractcristiquc dc 1' invention, on realise une correlation entre la 
densite minerale osseuse, DMO, obtenue d partir de cette image 2D et ledit 
parametre a evalue a partir de la meme image 2D suivant le modele 
mathematique : 

5 Cu' = bo + bi* exp (bz * DAAO) * a 

- avec bo, bi, bz qui sont des coefficients obtenus par regression non 
lineaire et la prediction de la contrainte ultime Cu de Tos. 
On determine une correlation entre le parametre a et un parametre 3D du 
reseau trabeculaire de Tos et un parametre 3D peut etre la densite de 
10 connect ivite Xv 

Pour supporter cette description, des dessins sont annexes et les figures qu'ils 
presentent, illustrent la description et montrent essentiellement des resultats 
obtenus. 

- figure 1 : courbe des resultats avec un modele lineaire, Cu en f onction de Cu* 

15 - figure 2 : courbe des resultats avec un modele exponentiel, Cu en f onction de 

- figure 3 : representation graphique de Cu en f onction de Xv 

- figure 4 : representation de f s en f onction de Xv 

- figure 5 :courbe de la f onction V(r) 

20 - figure 6 : determination du parametre a a partir de la courbe log-log de la 
f onction V(r) 

- figure 7 : courbe DMO en f onction de fs 

- figure 8 : courbe de a en f onction de Xv 

- figure 9 : courbe DMO en f onction de a 

25 - figure 10 : modele lineaire Cu en f onction de DMO 

- figure 11 : courbe de Cu en f onction Cu*. 

L' etude porte sur un os trabeculaire ou spongieux qui constitue environ 20% de 
la matiere osseuse. L'os cortical qui Tentoure assure la rigidite de I'os complet 
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ct done du squelette tandis que Tos trabcculaire assure la reprise et la 
resistance aux forces de compression. 

Selon {'invention, on considere que la structure trabeculaire comprend des 
travees qui sont interconnectees et la resistance mecanique de T ensemble 
5 rcsulte de cette conncctivite, c*est-a-dire du nombre de boucles fermees. La 
resistance mecanique resulte aussi du degre de mineralisation de ces travees, 
mineralisation qui est un parametre important. 

Le procede consiste a analyser une image numerisee obtenue dans le cas present, 
a partir d'une emission de rayons X. Cette image est une projection sur une 
10 surface d'une structure 3D ayant ete traversee par les photons emis. 

Chaque pixel qui const itue T image peut etre analyse independamment mais alors 
on ne peut determiner que la densite minerale osseuse DAAO si bien qu'il convient 
non seulement d 'analyser un pixel de fagon isolee mais de les analyser les uns par 
rapport aux autres. 

15 Le spectre du rayonnement doit etre le plus monochromatique possible af in 
d'eviter une trop forte derive dans Tenergie des photons. 

En eff et, si Tenergie varie la penetration variera et donnera des nuances de gris 
sur r image pouvant induire des erreurs en faisant croire a une absorption la ou il 
n'y en a pas. 

20 A partir de cette image obtenue avec un spectre etroit, on sait que chaque 
volume elementaire subit une action directe des photons, perpendiculairement 
mais aussi une action indirecte des rayons incidents diffuses par les volumes 
elementaires voisins apres qu'ils aient ete eux-memes traverses. 
Af'm de supprimer la part de diffusion, on peut utiliser de preference deux 

25 energies distinctes. On dispose de 2 equations a 2 inconnues et on peut done 
deduire la part d' attenuation reelle qui interesse directement le present 
procede et la part resultant de la diffusion. 



Si Ton recherche la contrainte ultime Cu (MPa), il faut tenir compte des 
parametres suivants : 

- Vo : volume d' analyse (mm^) 

- Vs : volume du solide (mm^) 

5 - f s ^ fraction volumique de solide 
avec la relation f g = Vs / Vo 

- Po • nombre d'amas de la partie solide (c*est-a-dire Tensemble des parties 
connectees ou disconnectees), ce nombre est generalement egal a 1. 

- P2 * nombre de surfaces internes, c'est-a-dire les trous resultant de la micro 
10 porosite interne des travees or d Techelle de travail, cette poroslte ne se voit 

pas, ce nombre est generalement egal a 0. 

- Nep nombre d'Euler-Poincarre 

- X : connectivite 

avec la relation x = Po + P2 - Nep 
15 - Xv • densite de connectivite (mm"^) 
avec la relation Xv = x / Vo 

Le procede selon la presente invention vise d permettre de correler deux 
parametres mesurables et de trouver la relation qui les lie. 
20 1 / Si Ton utilise un modele lineaire tel que : 

Cu* = Oo * fs + Cll 

- Cu' etant une valeur predite de Cu 

- oo et oi sont des coefficients de regression lineaire, 
il s'agit Id d'une simple correlation de Cuovec fs. 

25 La courbe obtenue est representee sur la figure 1 et Ton constate que la 
dispersion est tres importante. Si Ton se place par exemple d fs constant, 
I'ecart des valeurs de Cu est important. 
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On ne pcut rctenir ce modcle lincaire dont les rcsultats sont insuff isants en 
precision mais qui pcut ncanmoins servir dc comparatif . 

2/ Si Ton utilise un modele exponentiel tel que : 
5 Cu'= bo + bi* exp(b2 * fs) * Xv 

Avec bo, bi, bz qui sont des coefficients obtenus par regression non lineaire 
(methode des moindres carres generalisee) et Cu la prediction de la contrainte 
ultime Cu de Tos. 

La dispersion est fortement limitee, voir figure 2. Cu est estime a partir de fs 
10 etxv. 

Ces variations sont representees sur la figure 3. 

Si fs et Xv augmentent, Cu augmente. Ceci est parf aitement naturel, si la fraction 
volumique de solide augmente et la connectivite augmente, alors la contrainte 
ultime augmente. 

15 Si on fixe xv alors Cu est directement fonction de fs, de fagon exponentielle. A 
connectivite constante, la contrainte ultime augmente en fonction de 
r augmentation de la dens it e osseuse. 

Si on fixe fs alors Cu est directement fonction de Xv, de fagon lineaire mais la 
contrainte ultime Cu diminue lorsque la connectivite augmente. 
20 Si la fraction solide est constante, on constate que la resistance ultime Cu 
diminue lineairement en fonction de la connectivite. 

Plus on a de noeuds dans une structure et plus cette structure devient fragile, 
ceci de fa^on surprenante, contrairement a un prejuge bien etabli. 
On peut alors tracer le modele de representation de la figure 4 qui montre les 
25 valeurs de Cu. 

Pour la determination de la connectivite, il faut done trouver un parametre qui 
suive une meme loi et qui soit le plus independant possible de la densite osseuse 
pour ne pas etre influence. 
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On pourrait realiser un calcul a partir de deux parametrcs mcsurablcs comme la 
fraction volumique osseuse ct la conncctivitc, ccs deux paramctres ctant 
obtenus notammcnt par Imagerie a Resonance Magnetique. 

Par contre, Tobtention des valeurs par cette methode met en oeuvre des 
5 apparei'llages couteux, peu disponibles, ce qui conduit a des predictions d'un prix 
de revfent eleve interdisant un suivi longitudinal et done une diffusion massive 
aupres des patients. 

11 faut done pouvoir trouver deux valeurs mesurables par exemple a partir d'une 
simple densitometrie, qui est, elle, tout a fait abordable et qui peut etre 
10 renouvelee regulierement. En effet, des densitometries sont dejd realisees de 
fagon massive pour de nombreuses applications. La diffusion aupres des patients 
peut etre tres importante. 

Le procede selon la presente invention consiste a determiner ces deux 
parametres mesurables comme la fraction volumique osseuse et la connectivite, 

15 ou plus exactement des parametres qui pourraient leur etre correles de fagon 
suf f isamment etroite pour donner des resultats satisfaisants. 
Le premier parametre rapidement mesurable et f iable est la densite minerale 
osseuse. Des appareils ont ete developpes pour effectuer ces mesures de fagon 
f iable et reproductible. Si Ton trace la courbe DMO en fonction de fs, on 

20 constate que la dispersion est faible et que Ton peut utiliser DMO en lieu et 
place de f s.. Voir figure 7. Le premier parametre est obtenu. 
Le second parametre est appele a et determine a partir d'une image 2D 
numerisee par exemple obtenue avec un rayonnement X. 

Pour determiner ce nouveau parametre a, on choisit tout d'abord au hasard un 
25 pixel de T image 2D qui est a un niveau de gris h(0). 

On choisit ensuite une droite issue de ce point ayant une direction determinee 
egalement au hasard. 
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On se deplace d'une distance r sur ccttc droite, h(r) ctant Ic niveau de gris de 
ce nouveau point. 

On calcule alors la variance des niveaux de gris, ce qui donne : 

V(r)=[h(r)-h(0)]^ 

5 Le processus est iteratif , ovec une repetition de ces etapes un nombre de fois 
suff isamment grand, et aleatoire pour determiner a. On fait ainsi converger la 
fonction variance V(r) moyennee sur 1' ensemble des iterations. Ceci permet 
d'obtenir par calcul un ensemble de points et la courbe associee est representee 
sur la figure 5. 

10 Ainsi, plus on determine de points ayant un meme niveau de gris d proximite du 
pixel choisi, plus on se trouve dans la matiere et done plus on est homogene et 
plus on doit trouver de noeuds de connexion. 

Si Ton trace la courbe log-log de cette fonction, figure 6, on trouve sur les cinq 
premiers points en T occurrence, un comportement lineaire et la pente de cette 
15 droite est le coefficient denomme a et propre a la presente invention. 

Si Ton trace maintenant la courbe a en fonction de la connectivite volumique Xv, 
figure 8, on constate qu'il y a une correlation suffisante pour que Ton puisse 
remplacer Xv par a. 

On peut deduire le diagramme de la figure 9 de ces courbes precedentes, 
20 diagramme qui montre la relation entre DMO, a et Cu . 

DMO et a sont deux parametres directement mesurables et/ou pour le moins 
calculables a partir d'une meme image en 2D. 

On constate sur ce diagramme que lorsque la densite minerale osseuse reste 
constante, la contrainte ultime diminue lorsque a augmente. 
25 Ceci est bien equivalent au modele exponentiel precedemment indique avec la 
fraction osseuse en fonction de la connectivite, voir figure 4, sauf que la 
connectivite Xv n'etait pas mesurable directement alors que a est calcule a 
partir de T image 2D, done quantifiable. 
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On peut done ecnre la memc relation mais avec les parametrcs determinablcs a 
partir dc la meme image : 

Cu = bo + bi* exp(b2 * DMO) * a 
C'est la representation de la figure 9. La courbe de la figure 11 caracterise ce 
5 modele exponentiel de Cu' en fonction de a. On constate une faible dispersion 
comparee a celle du modele lineaire dont la courbe est mentionnee a la figure 10. 
Ceci est un exemple non limitatif qui montre la possibilite de lier la densite 
osseuse dont la valeur est determinee a partir d'une image 2D, a un parametre 
de structure determine a partir de cette meme image 2D pour obtenir la 
10 contrainte ultime Cu. 

On constate qu'il est possible de determiner une correlation entre le parametre 
a evalue a partir de T image 2D et une caracteristique 3D du reseau trabeculaire 
de Tos comme la densite de connectivite Xv. 

On peut ainsi donner un exemple obtenue a partir d' images osseuses femorales 
15 obtenue absorptiometrie a double energie qui remplace une image obtenue par 
rayons X ou &amma de faible energie. 

Ces images sont celles Tune d'une jeune femme de 17 ans et I'autre d'une 
femme agee de 69 ans. Les images 2D ainsi obtenues sont analysees suivant le 
precede de T invention. 
20 On constate a la mesure de la DMO les memes valeurs obtenues directement par 
le systeme densitometrique, dans les deux cas : 0,730 g/cm^ dans la zone 
trochanterielle. 

Les parametres a calcules sont respect ivement 0,56 pour le sujet jeune et 0,71 
pour le sujet age, montrant ainsi une difference notable et le parametre a le 
25 plus faible correspond a la meilleure resistance mecanique de I'os Cu, celui du 
sujet le moins age. 



12 

REVENDIC/\TIONS 

1. Procede de determination de la resistance mecanique d'un os a partir 
d*une image 2D numerisee, obtenue par imagerie, caracterise en ce que Ton 
realise une correlation entre la densite minerale osseuse DMO determinee a 
partir de cette image 2D par tout moyen adapte a ce type d'image et un 

5 parametre de structure obtenu a partir de la meme image 2D. 

2. Procede de determination de la resistance mecanique d'un os selon la 
revendication 1, caracterise en ce que Ton recourt a une correlation de type 
exponent ielle. 

3. Procede de determination selon la revendication 1 ou 2, caracterise en 
10 ce que la correlation associant la densite minerale osseuse et ledit parametre de 

structure est utilisee pour determiner la contrainte ultime Cy de I'os. 

4. Procede de determination de la resistance mecanique d'un os selon Tune 
quelconque des revendications precedentes, caracterise en ce que Ton determine 
le parametre de structure a obtenu par la succession des etapes suivqntes : 

15 a) choisir au hasard un pixel de 1* image 2D qui est a un niveau de gris h(0), 

b) choisir une droite issue de ce point ayant une direction egalement determinee 
au hasard, 

c) se deplacer d'une distance r sur cette droite, h(r) etant le niveau de gris de 
ce nouveau point, 

20 d) calculer la variance des niveaux de gris avec la formule : V(r) = [h(r) - h(0)]^, 

e) tracer la courbe associee a V(r) en echelle log-log, et 

f ) determiner la pente de cette courbe log-log qui represente ledit parametre a. 

5. Procede de determination de la resistance mecanique d'un os selon la 
revendication 4, caracterise en ce que Ton repete les etapes a) a d) un nombre 
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dc fois suffisamment grand pour fairc converger la fonction variance V(r), 
moyennee sur T ensemble des iterations. 

6. Procede de determination de la resistance mecanique d'un os selon la 
revendication 4 ou 5, caracterise en ce que Ton realise une correlation entre la 

5 densite minerale osseuse, DMO, obtenue a partir de cette image 2D et ledit 
parametre a evalue a partir de la meme image 2D suivant le modele 
mathematique : 

Cu' = bo + bi* exp (bz * DAAO) * a 
- avec bo, bi, bz qui sont des coefficients obtenus par regression non 
10 lineaire et Cu' la prediction de la contrainte ultime Cu de I'os. 

7. Procede de determination de la resistance mecanique d'un os selon Tune 
quelconque des revendications 4, 5 ou 6, caracterise en ce que Ton determine 
une correlation entre le parametre a et un parametre 3D du reseau trabeculaire 
de Tos. 

15 8. Procede de determination de la resis-tance mecanique d'un os selon la 

revendication 7, caracterise en ce que le parametre 3D du reseau trabeculaire 
de I'os est la densite de connect ivite Xv 
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ABREGE DESCRIPTTF 



L'objet de 1' invention est un procede de determination de la resistance 
mecanique d'un os a partir d'une image 2D numerisee, obtenue par imagerie, 
caracterise en ce que Ton realise une correlation entre la densite minerale 
osseuse DAAO determinee a partir de cette image 2D par tout moyen adapte a ce 
type d'image et un parametre de structure obtenu a partir de la meme image 2D. 



